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ABSTRACT 

We report the first Chandra detection of emission out to the virial radius in the 
cluster Abell 1835 at z = 0.253. Our analysis of the soft X-ray surface brightness 
shows that emission is present out to a radial distance of 10 arcmin or 2.4 Mpc, and 
the temperature profile has a factor of ten drop from the peak temperature of 10 keV 
to the value at the virial radius. We model the Chandra data from the core to the virial 
radius and show that the steep temperature profile is not compatible with hydrostatic 
equilibrium of the hot gas, and that the gas is convectively unstable at the outskirts. 
A possible interpretation of the Chandra data is the presence of a second phase of 
warm-hot gas near the cluster's virial radius that is not in hydrostatic equilibrium 
with the cluster's potential. The observations are also consistent with an alternative 
scenario in which the gas is significantly clumped at large radii. 

Key words: galaxies: clusters: individual (Abell 1835); cosmology: large-scale struc- 
ture of universe. 



1 INTRODUCTION 

The large-scale halo of hot gas provides a unique way to 
measure the baryonic and gravitational mass of galaxy clus- 
ters. The baryonic mass can be measured directly from 
the observation of the hot X-ray emitting intra-cluster 
medi um (ICM), and of the associated stella r component 
(e.g. iGiodini et all 120091 ; iGonzalez et all 120071 ') . while mea- 
surements of the gravitational mass require the assumption 
of hydrostatic equilibrium between the gas and dark matter. 
Cluster cores are subject to a variety of non-gravitational 
heating and cooling processes that may result in deviations 
from hydrostatic equilibrium, and in inner regions beyond 
the core the ICM is expected to be in hydrostatic equilibrium 
with the dark matter potential. At the outskirts, the low- 
density ICM and the proximity to the sources of accretion re- 
sults in the onset of new physical pr ocesses such as d eparture 
from hydrost atic equilibrium (e.g.. | Lau et al.ll2009h . clump- 



ing of the gas |Simionescu et al.l2011 ), different temper ature 
between electrons and ions (e.g..lAkamatsu et al. Il201ll ). and 
flattening of the entropy profile (|Sato et al.l 120121 ). leading 
to possible sources of systematic uncertainties in the mea- 
surement of masses. 

The detection of hot gas at large radii is limited 
primarily by its intrinsic low surface brightness, uncer- 
tainties associated with the subtraction of background 
(and foreground) emission, and the ability to remove 



contamination from compact sources unrelated to the 
cluster. Thanks to its low detector background, Suzaku 
reported the measurement of ICM temp eratures to r2oo 
and b e yond for a few nearby clus ters (e.g. lAkamatsu et al.l 



| 201ll; IWalker et al.l l2012bllal; ISimionescu et all 



2011 



burns et all 120101; I Kaw aharada et al.l |2010| ; IfBautz et al.l 



120091 ; iGeorge et ail 120091 ); to date Abell 1835 has not been 
the target of a Suzaku observation. 

In this paper we report the Chandra detection of X-ray 
emission in Abell 1835 beyond T2oo, using three observa- 
tions for a total of 193 ksec exposure time, extending the 
analys is of these Chandra data performed by I Sanders et al.l 
(2010). The radius ta is defined as the radius within which 
the average mass density is A times the critical density of 
the universe at the cluster's redshift for our choice of cos- 
mological parameters. The virial radius of a cluster is de- 
fined as the equilibrium radius of the collapsed halo, ap- 
proximat ely equivalent to one half of its tu rnaround ra- 
dius (e.g. lLacev fc Colelll993l ; lEke et al.lll99Sl ). For an n A - 
domi nated universe, the virial radius is approximately rioo 
(e.g. lEke et al.ll 199Sf) . A bell 1835 is the most luminous clus- 
ter in the lDahlel (|2006h sample of clusters at z = 0.15 — 0.3 
selected from the Bright Cluster Survey. The combina- 
tion of high luminosity and availability of deep Chandra 
observations with local background make Abell 1835 and 
ideal candidate to study its emission to the virial radius. 
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Abell 1835 has a redshift of z = 0.253, which for H = 
70.2 km s" 1 Mnc" 1 ft A = 0.73, Cl M = 0.27 cosmology 
l|Komatsu et al.l l201ll ) corresponds to an angular-size dis- 
tance of Da = 816.3 Mpc, and a scale of 237.48 kpc per 



2 CHANDRA AND ROSAT OBSERVATIONS 
OF ABELL 1835 AND THE DETECTION OF 
CLUSTER EMISSION BEYOND r 20 o 

2.1 Chandra observations 

Chandra observed Abell 1835 three times between Decem- 
ber 2005 and August 2006 (observations ID 6880, 6881 and 
7370), with a combined clean exposure time of 193 ks. The 
three observations had similar aimpoint towards the cen- 
ter of the cluster (R.A. 14h01m02s, Dec. +02d51.5m J2000) 
and different roll angles. All observations were taken with 
the ACIS-I detector configuration, which consists of four 
ACIS front-illuminated chips in a two-by-two square, plus 
a fifth identical chip that may be used to measure the 
in situ soft X-ray background. Figure [1] is an image from 
the longest observation (ID 6880, 118ks) in the soft X- 
ray band (0.7-2 keV). In addition to a large number of 
compact X-ray sources that were excluded from further 
analysis, the data show a clear detection of diffuse X-ray 
emission associated with two additional low-mass clusters 
identified from the Sloan Digital Sky Survey, MAXBCG 
J210. 31728+02. 75364 and WHL J140031.8+025443. The 
cluster MAXBCG J210. 31728+02. 75364 is the only cluster 
in th e vicinity of Ab e ll 183 5 reported in the MAXBCG cata- 
log of lKoester et al.l l|2007h, and it has a me asured photo- z of 
0.238, while the catalog of lWen et alj (|2009h reports a photo- 
z of 0.269 for the same source; given the uncertainties asso- 
ciated with photometric redshifts, it is likely that the clus- 
ter is in physical ass ociation with Abell 1835 (z = 0.253). 
The lWen et ail l|2009h catalog also reports another optically- 
identified cluster in the area, WHL J140031.8+025443, 
with a spectroscopic redshift of z = 0.2505. The associ- 
ation of these two groups with Abell 1835 is confirmed 
by redshift data provided by C. Haines (personal com- 
munication), who measures a redshift of z = 0.250 for 
WHL J140031.8+025443, and z = 0.245 for MAXBCG 
J210.31728+02.75364. 

Since the goal of this paper is to study the diffuse emis- 
sion associated with Abell 1835, we excise a region of radius 
90 arcsec around the position of the two clusters (black cir- 
cles in Figure [TJ, and study their emission separately from 
that of Abell 1835 (see Section E2J . 



2.2 Chandra data analysis 



C handra observation s follo w the 
Bonamente et all ([2006 ) and 



The reduction of the 
procedure described 

iBonamente fc Nevalainen! l |201ll ). which consists of filtering 
the observations for possible periods of flaring background, 
and applying the latest calibration; no significant flares were 
present in these observations. The reduction was performed 
in CIAO 4.2, using CALDB 4.3; in Sec. E3]we discuss the 
impact of calibration changes on our results. One of the cal- 
ibration issues that can affect the measurement of cluster 
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Figure 1. Image of Abell 1835 from observation 6880, in the 0.7- 
2 keV band. The data were smoothed with a Gaussian kernel of 
~6 arcsec standard error. The dashed circles correspond to radial 
distances of approximately rsoo and r200, and the full black circles 
mark the position of the two low-mass clusters associated with 
Abell 1835. 



emission is the uncertainty in the contamination of the opti- 
cal blocking filter, which causes a reduction in the low energy 
quantum efficiency of the Chandra detectors. The spatial 
and time dependence of this contaminant affects primarily 
the effective area at ^0.7 keVQ with an estimated residual 
error of ^ 3% at higher energy. We therefore limit our spa- 
tial and spectral analysis to the ^0.7 keV ba nd. The superior 
angu lar resolution of the Chandra mirrors dWeisskopf et al.l 
l2000h results in a point-spread function with a 0.5 arcsec 
FWHM, and therefore there is negligible contribution from 
the bright cluster core to the emission in the outer annuli, 
and from secondary scatter (stray light) by sources outside 
the field of view. 

The subtraction of particle and sky background is one 
of the most crucial aspects of the analysis of low surface 
brightness cluster regions. We use Chandra blank-sky back- 
ground observations, rescaled according to the high- energy 
flux of the cluster, to ensure a correct subtraction of the par- 
ticle background that is dominant at E ^ 9.5 keV, where the 
Chandra detectors have no effective area. The temporal and 
spatial variability of the soft X-ray background at E < 2 keV 
also requires that a peripheral region free of cluster emission 
is used to measure any local enhancement (or deficit) of soft 
X-ray emission relative to that of the blank-sky fields, and 
account for this difference in the analysis. This method is ac- 
curate for the determination of the temperature profile, but 
may result in small errors in the measurement of the surface 
brightness profile. In fact, the blank-sky background is a 
combination of a particle component that is not vignetted, 
and a sky component that is vignetted. To determine the 
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Figure 2. Exposure corrected surface brightness profile of 
Abell 1835 in the soft X-ray band (0.7-2 keV), obtained by sub- 
traction of the particle background from the ACIS stowed obser- 
vations. The radii rsoo, »~200 and the virial radius (~ rioo) are 
estimated from the data in Section l4.2l (see Table |4j. The dashed 
red line is the average background level in the region ^ 700 arcsec. 

surface brightness of the cluster and of the local soft X- 
ray background, a more accurate procedure consists of sub- 
tracting the non-vignetted particle component as measured 
from Chand ra observations in which the ACIS detector was 
stowed (e.g.. lHickox fc Markevitchll2007n . after rescaling the 
stowed background to match the E ^ 9.5 keV cluster count 
rate, as in the case of the blank-sky background. 

Point sources are identified and removed using a wavelet 
detection method that correlates the cluster observation 
with wavelet functions of different scale sizes (wavdetect in 
CIAO). Subtraction of point soures from the blank-sky ob- 
servations were performed by eye, with results that closely 
match those of the wavelet method. 

2.3 Measurement of the surface brightness profile 
with Chandra 

The surface brightness profile obtained using this back- 
ground subtraction is shown in Figure [2] in which the red 
line represents the average value of the background at radii 
700 arcsec, where the surface brightness profile is consis- 
tent with a constant level. To determine the outer radius at 
which Chandra has a significant detection of the cluster, we 
also include sources of systematic errors in our analysis. One 
source of uncertainty is the error in the measurement of the 
background level, shown in Figure [3] as the solid red lines. 
The error is given by the standard deviation of the weighted 
mean of the datapoints at radii greater than 700 arcsec, to 
illustrate that each bin in the surface brightness profile be- 
yond this radius is consistent with a constant level of the 
background. 

Another source of uncertainty is the amount by which 
the stowed background is to be rescaled to match the cluster 
count rate at high energy. The stowed background dataset 
applicable to the dates of observation of Abell 1835 has 
an exposure time of 367 ksec, and the relative error in 
the rescaling of the background to match the cluster count 
rate at high energy is 0.7%, as determined by the Pois- 
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Figure 3. Close-up view of Figure[2] in which the red lines rep- 
resent the l-cr confidence in the background level as determined 
from the ^700 arcsec region, and the green error bars combine 
the statistical and systematic errors in the determination of the 
surface brightness. 

son error in the photon cou nts at high energy. Moreover, 
iHickox fc Markevitchl (2006) has shown that the spectral 
distribution of the particle background is remarkably sta- 
ble, even in the presence of changes in the overall flux, and 
that the ratio of soft-to-hard (2-7 keV to 9.5-12 keV) count 
rates remains constant to within ^2 %. We therefore apply 
a systematic error of 2 % in the stowed background flux, to 
account for this possible source of uncertainty, in addition 
to the 0.7% error due to the uncertainty in the rescaling of 
the background. 

In Figure [3] the green error bars represent the cumula- 
tive effect of the statistical error due to the counting statis- 
tic, and the sources of errors associated with the use of the 
stowed background; the systematic errors were added lin- 
early to the statistical error as a conservative measure. This 
error analysis shows that the emission from Abell 1835 re- 
mains significantly above the background beyond r2oo and 
until approximately a radius of 600 arcsec, or approximately 
2.4 Mpc. The significance of the detection in the region 450- 
600" (the five datapoints in Figure [3] after the r2oo marker) 
is calculated as 5.5<r, and is obtained by using the larger sys- 
tematic error bars for the surface brightness profile (in green 
in Figure [3]), added in quadrature to the error in the deter- 
mination of the background level from the 700" region 
(red lines in Figure [3]). 

To further test the effect of the background subtrac- 
tion, we repeat our backround subtraction process using the 
^ 600" region (instead of the 700" region) . The back- 
ground level increases by less than la of the value previ- 
ously determined (e.g., the two levels are statistically indis- 
tinguishable), and the significance of detection in the region 
450-600" is 4.7cr. Therefore we conclude that it is unlikely 
that the excess of emission beyond r2oo and out to the virial 
radius is due to errors in the background subtraction pro- 
cess. A similar result can be obtained including the 2-7 keV 
band, but the signal-to-noise is reduced because at large 
radii this band is dominated by the background due to the 
softening of the cluster emission. We estimate r2oo and the 
virial radius (~ rioo) from the Chandra data in Section \4. 21 
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Figure 4. Surface brightness profile in 0.2-2 keV band from a 
6 ks observation with ROSAT PSPC. The background level is 
determined from the data at radii 700", as in the Chandra 
data. 




0.5 1 2 5 



Energy (keV) 

Figure 5. Spectrum of the local enhancement of the soft X-ray 
background from observation 6880. The other two exposures have 
similar levels of soft X-ray fluxes above the blank-sky emission, 
which is modeled as an unabsorbed ~ 0.25 keV thermal plasma 
at z = 0. The best-fit model has a Xmin = 73.9 for 78 degrees of 
freedom, for a null hypothesis probability of 61%. 



2.4 Measurement of the surface brightness profile 
with the ROSAT Position Sensitive 
Proportional Counter 

R OSA T observed Abell 1835 on July 3-4 2003 for 6 ks with 
the Position Sensitive Proportional Counter (PSPC), obser- 
vation ID was 800569. The PSPC has a 99 .9% rejection of 
parti cle background in the 0.2-2 keV band (|Plucinskv et al.l 
1 19931 ) and an average angular resolution of ~30 arcsec that 
makes it very suitable for observations of low surface bright- 
ness objects such as the outskirt s of galaxy clusters (e.g. 
IBonamente et al.ll200ll , |2002| . I2003D . We reduce the event file 
follo wing the procedu r e desc ribed in ISnowden et all (| 19941 ) 
and IBonamente et all l|2002l ) , which consists of corrections 
for detector gain fluctuations, and removal of periods with 
a master veto rate of ^170 counts s in order to discard 
periods of high background. These filters result in a clean 
exposure time of 5.9 ks. 

Since the PSPC background is given only by the photon 
background, we generate an image in the 0.2-2 keV band and 
use the exposure map to correct for the position-dependent 
variations in the detector response and mirror vignetting. 
We masked out the two low-mass cluster regions as we did 
for the Chandra data and all visible point sources, and ob- 
tained and exposure-corrected surface brightness profile out 
to a radial distance of ~20 arcmin, which corresponds to the 
location of the inner support structure of the PSPC detec- 
tor. The ROSAT surface brightness profile therefore covers 
the entire azimuthal range. In Figure]?] we show the radial 
profile of the surface brightness in the 0.2-2 keV band, show- 
ing a ~2 a excess of emission in the 400-600" region using 
the background level calculated from the region ^700", as 
done for the Chandra data. The ROSAT data therefore pro- 
vide additional evidence of emission beyond rsoo and out to 
the virial radius, although the short ROSAT exposure does 
not have sufficient number of counts to provide a detection 
with the same significance as in the Chandra data. 



3 ANALYSIS OF THE CHANDRA SPECTRA 

3.1 Measurement of the temperature profile of 
Abell 1835 

We measure the temperature profile of Abell 1835 follow- 
ing the background subtraction method described in Sec. [2] 
which makes use of the blank-sky background dataset and 
a measurement of the local enhancement of the soft X-ray 
background, as is commonly done for Chandra data (e.g. 
Vikhlinin etall 120061 : iMaughan et all 120081 ; iBulbul et all 
20101 ). In Figure [51 we show the spectral distribution of the 
local soft X-ray background enhancement, as determined 
from a region beyond the virial radius (^700 arcsec); this 
emission was modelled with an APEC emission model of 
kT ~ 0.25 KeV and of Solar abundance, consistent with 
Galactic emission, and then subtracted from all spectra. The 
spectra were fit in the 0.7-7 keV band using the minimum 
X 2 statistic, after binning to ensure that there are at least 
25 counts per bin. We use XSPEC version 12.6.0s for the 
spectral analysys. 

In Figure we show the spectra of the outermost two 
regions, to show the impact of the soft X-ray residuals in 
the background subtraction. The importance of background 
systematics in the detection of emission and measurement 
of cluster temperatures fo r regions of low surfa c e brig htness 
was recently addressed bv lLeccardi fc Molendil (120081 ) using 
XMM-Newton data. For our Chandra observations, the two 
main sources of uncertainty when determining the tempera- 
ture of the outer regions are the subtraction of the blank-sky 
background, and the subtraction of the locally-determined 
soft X-ray background. Table [1] reports the statistics of the 
background relevant to the outer regions of the cluster, with 
both regions ~ 10-20% above the blank-sky background, de- 
termined with a precision of 1-2%. The additional soft X-ray 
background accounts for a significant portion of the remain- 
ing signal, as shown in Figure HJ] the 90% upper limit to the 
measurement of this background is shown as the green lines, 
and emission from the cluster is still detected with high sta- 



te) 2012 RAS, MNRAS 000, [11,131 



Chandra X-ray observations of Abell 1835 to the virial radius 5 



o - 
b 



> 




Energy (keV) Energy (keV) 

Figure 6. Blank-sky background subtracted spectra of regions 330-450" and 450-600" from observation 6880. The solid lines are the 
best-fit model of the local soft X-ray enhancement of Figure \E\ (red), and its 90% upper limit (green). 



Table 1. Background levels in outer regions 



Exposure time (ks) 

Correction to Blank-sky Subtraction" 



Observation ID 
6880 6881 7370 



114.1 36.0 39.5 

-0.04±0.01 -0.125±0.015 -0.04±0.015 











Re f 


;ion 330-450" 






Total Counts 










18,124 


4,938 


5,686 


Count rate (c s — 1 ) 










0.158±0.001 


0.137±0.002 


0.144±0.002 


Net count rate 6 (10" 


-2 


c s 






2.75 ± 0.10 


2.54 ±0.20 


2.36 ±0.20 


Percent above back. 










17.4±0.6 


18.2±1.4 


16.4±1.4 


SXB count rate (10" 


3 


c s" 


x ) 




3.34±0.77 


7.20±0.94 


7.10±0.71 










Re; 


;ion 450-600" 






Total Counts 










15,811 


4,901 


5,483 


Count rate (c s — 1 ) 










0.139±0.001 


0.136±0.002 


0.139±0.002 


Net count rate 6 (10" 


-2 


c s~ 


x ) 




1.02±0.10 


1.58±0.20 


1.23±0.20 


Percent above back. 










7.3±0.7 


11.6±1.5 


8.8±1.4 


SXB count rate (10" 


-3 


c s~ 






3.06±0.70 


7.50±0.98 


7.32±0.73 



a: This is the fractional correction of the blank-sky data, to match the high-energy flux in the cluster observation. 
b This is the background-subtracted count rate, including cluster and soft X-ray background (SXB) signal. 



tistical significance. Both sources of error are included in the 
temperature measurements at large radii. 

We use the APEC code l|Smith et alj|200ll . code version 
1.3.1) to model the Chandra spectra, with a fixed Galactic 
HI co lumn density of Nh = 2.04 x 10 20 cm" 2 (|Kalberla et all 
120051 ). The region at radii 330" have a variable metal 
abundance, while the outer region have a fixed abundance 
of A = 0.3. In addition to the statistical errors obtained 
from the XSPEC fits, we add a systematic error of 10% in 
the temperature measured in the core and a 5% error to the 
other region, to account for possible systemat ic uncertain- 
ties d ue to the Chandra calibration (see, e.g., iBulbul et alj 
I2010T) . One possible source of systematic uncertainty in our 
results is indicated by the systematic difference between the 
Chandra/ ACIS andXMM-TVe wfon/EPIC temperature mea- 
surements of galaxy clusters (Nevalain en et ai1l2010l ) This 
amounts to a ±10% bias in the calibration of the effective 
area at 0.5 keV, which decreases roughly linearly towards 0% 



bias at 2 keV. Assuming that XMM-Newton/pn has a more 
accurately calibrated effective area, we reduced the Chan- 
dra effective area by multiplying it with a linear function as 
indicated by the Chandra/ XMM-Newton comparison. As a 
result, the temperature at the outermost radial bin decreases 
by ~ 5%. Thus, the cross-calibration uncertainties between 
Chandra and XMM-Newton do not explain the low tempera- 
ture we measure in the outermost radial bin. Uncertainties in 
the Galactic column density of HI do not impact significantly 
our results. Changing the value of Nh by ±10% , consi stent 
with the variations between the iKalberla et ail l|2005l ) and 
the iDickev fc Lockmanl |l990T l measurements, results in a 
change of best-fit temperature in each bin by less than 2%. 

Given the emphasis of this paper on the detection of 
emission at large radii, we investigate the sources of uncer- 
tainty caused by the background subtraction in the outer 
region at ^330". We report the results of this error analy- 
sis in Table [2l where cornorm refers to the normalization of 
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Table 2. Temperature measurement and error analysis from the 
Chandra data. 



Region Projected Temperature (keV) 





Measurement 11 


Calibration error' 


0-10" 


4.78±0.06 


±0.48 


10-20" 


7.09±0.14 


±0.71 


20-30" 


8.72±0.27 


±0.87 


30-60" 


9.47±0.21 


±0.47 


60-90" 


10.57±0.33 


±0.53 


90-120" 


9.97±0.44 


±0.50 


120-180" 


9.68±0.49 


±0.48 


180-240" 


7.85±0.65 


±0.39 


240-330" 


6.02±0.65 


±0.30 


330-450" 


3.75±0.72 


±0.19 


450-600" 


1.26±0.16 


±0.06 


Measurement of Temperature Using 


Background Systematic Errors (keV) 




+ lcr cornorm c 


— lcr cornorm 


330-450" 


3.02±0.54 


4.67±1.00 


450-600" 


1.09±0.10 


1.31±0.18 




+lcr soft res. d 


— la soft res. 


450-600" 


4.53±1.03 


3.05±0.54 


450-600" 


1.37±0.25 


1.16±0.12 


Summary of Background Systematic Errors 6 


330-450" 


±0.83 ±0.74 keV 


450-600" 


±0.11 ±0.10 keV 



a: Uncertainty is la statistical error from counting statistics 
only. 

b: Includes XMM-N ewton/ Chandra cross-ca libration uncertainty 

of the effective area {Nevalainen et alJlipToT) . 

c: This is temperature obtained by varying by ±lc the 

fractional correction of the blank-sky data, to match the 

high-energy flux in the cluster observations. 

d: This is the temperature obtained by varying by ±1<t the 

normalization of the best-fit model to the soft X-ray background 

residuals. 

e: Obtained from the average deviation of the ±1<t 'cornorm' 
and 'soft, res' measurements from the measurement with 
nominal values of these parameters. 

the blank-sky background, and soft residuals refers to the 
normalization of the soft X-ray residual model, as reported 
in Table [1] In the analysis that follows, we add the system- 
atic errors caused by these sources linearly to the statistical 
error. Our data do not constrain well the metal abundance 
of the plasma in the outer regions. Using an abundance of 
A = 0.5 instead of the nominal A = 0.3 leads to negligible 
changes in the best-fit temperature for both of the outer an- 
nuli. In the extreme case of an A — 0.0 metal abundance, 
both regions have an acceptable fit with the best-fit tem- 
peratures change respectively by +6% for the 330-450" re- 
gion (Ax 2 = +1.3), and by -22% for the 450-600" region 
(Ax 2 = +9.2, best fit decreases from 1.26 to 0.98 keV). We 
therefore find that, in the case of exceptionally low metalllic- 
ity, the temperature profile we measure from these Chandra 
data would be even significantly steeper than indicated by 
the result in Table [5] Given that these data do not provide 
direct indication that the plasma in the outer regions may 
have null metal content, we do not fold in this source of 
syst ematic error in the an alysis that follows. 

ISanders et ail l|2010l ) measured temperature profiles 
for Abell 1835 out to approximately rsoo with Chan- 



dra and XMM-Newton. Using the same Chandra obser- 
vations we analyze in this paper, their temperature pro- 
file has a similar drop from the peak value to their out- 
ermost annulus (322 ± 42"), where they measure a tem- 
perature of kT=4.67±o;H keV that is consistent with 
our measurements. Likewise, from the XMM-Newton data 
their outermost radial bin (300 ± 10") has a tempera- 
ture of kT=5.2±J; 2 keV, also in agreement with our re- 
sults. The only measurement of the Abell 1835 tempera- 
ture to the viri a l radi us available in the literature is that of 
ISnowden et~ai1 (|2008l ), who does report a temperature pro- 
file out to a distance of 12 arcmin from a long XMM-Newton 
observation (and out to 7' from a shorter observation). In 
particular, they report a temperature of kT — 3.14 + 0.93 
for the region 420-540", which straddles our measurements 
at 330-450" (3.75 ± 0.72 keV) and at 450-600" (1.26 ± 0.16, 
statistical errors only). The same paper also reports a mea- 
surement of kT = 3.33 ± 1.75 keV for the region 540-720", 
i.e., beyond our outer annulus. Their temperature is some- 
what higher that ours, although the large error bars cannot 
exclude that the Chandra and XMM-Newton measurements 
are consistent. Therefore our res ults confirm and exten d the 
earlier XMM-Newton analysis of ISnowden et all l|2008h . 

3.2 Measurement of the average temperature of 
MAXBCG J210. 31728+02. 75364 and WHL 
J140031. 8+025443 

We also measure the temperature of the two SDSS 
clusters detected in our Chandra images, MAXBCG 
J210. 31728+02. 75364 and WHL J140031.8+025443. The 
two clusters are located between a distance of ~380-650" 
from the cluster center, and therefore we start by extracting 
a spectrum for this annulus excluding two regions of 1.5' 
radius centered at the two clusters. This radius was deter- 
mined by visual inspection, after smoothing of the Chandra 
image with a Gaussian kernel of a = 6 arcsec. For this annu- 
lus, we measure a temperature of kT — 1.85 ± 0.36 keV for 
a fixed abundance of A = 0.3 Solar. We then use this spec- 
trum as the local background for the two cluster regions, and 
measure a temperature kT = 2.73±8:§4 keV i for MAXBCG 
J210. 31728+02. 75364 (357 source photons, 19% above the 
average emission of the annulus), and kT = 2.09+o:| 5 keV 
for WHL J140031.8+025443 (538 photons, 27% above back- 
ground). For both clusters, we assumed the same Galactic 
HI column density as for Abell 1835, and a fixed metal abun- 
dance of A = 0.3 Solar. For both clusters we also extract 
spectra in regions larger than 1.5', and determine that no 
additional source photons are present from these two clus- 
ters beyond this radius. 

3.3 Tests of robustness of the temperature 
measurement at large radii 

To further test the measurement of temperatures especially 
at large radii, where the background subtraction is espe- 
cially important, we also measure the temperature profile 
using the same stowed background data that was used for 
the surface brightness measurement of Figures [2] and [3] As 
in the case of the blank-sky background, we first rescale the 
stowed data to match the high-energy count rate of the clus- 
ter observation, and use a region at large radii (^ 700 arcsec) 
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to measure the local X-ray background. We model the back- 
ground using an APEC plus a power-law model, the latter 
component necessary to model the harder emission due to 
unresolved AGNs that is typically removed when the blank- 
sky background is used instead, and apply this model to all 
cluster regions. We find that the temperature profile is con- 
sistent within the 1 a statistical errors of the values provided 
in Table [2] for each region, and therefore conclude that the 
temperature drop at large radii, and especially in the outer- 
most region, is not sensitive to the background subtraction 
method. 



The temperature measurement is also dependent on 
an accurate subtraction of background (and foreground) 
sources of emission. Point sources in the field of view are de- 
tetected using the CIAO tool wavdetect, which correlates the 
image with wavelets at small angular scales (2 and 4 pixels, 
one pixel is 1.96"), searches the results for 3-a correlations, 
and returns a list of elliptical regions to be excluded from 
the analysis. We study in particular the effect of background 
sources on the measurement of the temperature in the out- 
ermost annulus (450-600"). In this region, wavdetect finds 
24 point sources, plus portions of the two low-mass galaxy 
clusters described in Section T3. 21 We extract a spectrum for 
this region from the longest observation (ID 6880) , and now 
include in the spectrum all point sources excluded in the 
previous analysis. We find a count rate of 3.20 ± 0.13 x 10 -2 
counts s _1 , compared to the point source-subtracted rate 
of 1.02 ± 0.10 x 10~ 2 counts s~\ corresponding to an in- 
crease in background-subtracted flux by a factor of three. 
We then fit the spectrum with the same APEC model as 
described in Section [3.11 and find a best-fit temperature of 
kT — 1.96 ± 0.17 keV for a best-fit goodness statistic of 
X 2 = 537 for 429 degrees of freedom (or x^. ed =1.25), com- 
pared to the temperature of 1.22 ± 0.19 keV for a \ 2 = 415 
for 389 degrees of freedom (or x 2 e d=l-08). We therefore con- 
clude that an accurate subtraction of point sources and un- 
related sources of diffuse emission is crucial to obtain an 
accurate measurement of the temperature profile, especially 
in regions of low-surface brightness such as those near the 
virial radius. 



Changes in the instrument calibration affect the mea- 
surement of temperatures. We therefore repeat the same 
data reduction and spectral analysis using the latest soft- 
ware and calibration database available at time of writing 
(CIAO 4.4 and CALDB 4.5.1) for the longest observation 
(ID 6880) , and obtain a new temperature profile for the same 
regions as reported in Table [2] In the outermost two regions, 
we measure a temperature of 3.04±0.69 keV (330-450") and 
1.23±0.21 keV (450-600"), well within the 1-a confidence 
intervals of the measurements using the older calibration 
(3.40±0.76 and 1.22±0.19 respectively, also in agreement 
with the values of Table [2] obtained from the combination 
of all exposures). The temperature of the inner regions are 
also always within 1-a of the results obtained with the earlier 
calibration, and we therefore conclude that changes in the 
instrument calibration do not affect significantly our results. 



4 MEASUREMENT OF MASSES AND GAS 
MASS FRACTION 

We fit t he surface brigh t ness and the temperature profiles 
with the lVikhlinin et all (|2006l ) model. The electron density 
is modelled with a double-/? profile modified by a cuspy core 
component and an exponential cutoff at large radii, for a to- 
tal of eleven model parameters; the temperature has both a 
cool-core component to follow the cooler gas in the core, and 
a decreasing profile at large radii, for an additional n i ne pa - 
rameters. For our analysis, we follow Ivikhlinin et all (2006) 
and fix the 7=3.0 parameter, and do not use the cuspy-core 
component (a = 0) or the second /3-model component, so 
that the density is modelled by just one /3-model with an 
exponential cutoff, for just four free parameters (core radius 
r c , exponent /3, scale radius r a and exponential cutoff expo- 
nent e, see Table [3}. For the temperature profile, we fix the 
parameter a = 0, and the remaining eight parameters are 
reported in Table [3] 

We use a Monte Carlo Markov chain (MCMC) method 
that we us ed in previous papers (e.g., iBonamente et all 
12004 120061 ). The MCMC analysis consists of a projection 
of the three-dimensional models and a comparison of the 
projected surface brightness and temperature profiles, and 
results in simultaneous estimation of the posterior distribu- 
tions of all model paramters. Uncertainties in the parameters 
are obtained from the posterior distributions, with 1-a er- 
rors assigned using the 68.3% confidence interval around the 
median of the distribution. 

The gas mass is directly calculated from the electron 
density model parameters via 



M ga3 (r) = m p /j, e / n e (r)4nr dr 



(1) 



and the total gravitational mass via the equation of hydro- 
static equilibrium, 



kT(r)r f d\nn e dlnkT 

Mir) = H 

/j, e m p G \ dlnr dlnr 



(2) 



where m p is the proton mass, [i e ~ 1.17 the mean electron 
molecular weight, and G the gravitational constant. The to- 
tal density of matter is simply obtained via 



p(r) 



1 dM(r) 



47T7" 2 dr 

and therefore can be obtained via a derivative of the 
mass profile. In Equation [2] the term A = din n e / dlnr + 
dlnkT /dlnr and its first derivative are always negative, as 
is dkT(r)/dr at large radii. Therefore, the density can be 
rewritten as 



p(r) = 



47rr 2 fi e m p G 



kT 



dA 

A + r — 

dr 



+ rA 



dkT 
dr 



(3) 



in which the only negative term is the one containing A ■ 
dkT(r) /dr, while the other two terms remain positive out 
to large radii. 



4.1 Modelling of the Chandra data out to the 
virial radius 



The lVikhlinin et al.l (|2006l ) model provides a satisfactory fit 
out to the outermost radius of 600" ; Figure [7] shows the 
best-fit models to the temperature and surface brightness 
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profiles, best-fit parameters of the model are reported in 
Table [3] The temperature profile measured by Chandra in 
Figure [7] is so steep that it causes the total matter density 
p(r) to become negative at approximately 400", indicating 
that the temperature profile cannot originate from gas in 
hydrostatic equilibrium. The situation is illustrated in Fig- 
ure where the relevant terms of Equation [3] are plotted 
individually; the density inferred from hydrostatic equilib- 
rium becomes negative where the negative term crosses the 
positive ones, and the mass profile has a negative slope be- 
yond that point. These fit parameters therefore lead to an 
unacceptable situation, and responsibility for this inconsis- 
tency can be attributed to an overly steep temperature pro- 
file, with a drop by a factor of ten between approximately 
1.5' to 10'. 

The results presented in this section provide evidence 
that the gas detected by Chandra near the virial radius is 
not in hydrostatic equilibrium, and a number of theoreti- 
cal studies do in fact suggest that beyond rsoo the inter- 
galactic p lasma is not su pported solely by thermal pres- 
sure (e.g. lLau et all 120091 ') . Suzaku has reported the mea- 
surement of emission near the virial radius for several clus- 
ters, including Abell 1413, Hydra A, Pe rseus, PKS0745-191 
Abell 1795, Abell 1689 and Abell 2029 dHoshino et al.ll201C 



Sato et al.ll2012l;ISim ioncscu et al. 2011 ; iGeorge et al.l 2009; 
Bautz et all 120091; iKawaharada et all l20ld; Walker et all 



20123E). Some of these results do in fa ct report an appar- 



ent decrease in total ma ss with radius (George et al.l 12009 ; 
IKawaharada et aT]|2010l . e .g.) and lack of hy drostatic equi- 
librium at large radii fe.g. iBautz et al.1l2009l ). similar to the 
results presented in this paper. Temperature profiles mea- 
sured by Suzaku typically do not feature as extreme a tem- 
perature drop as the one reported in Figure [7] i.e., a factor 
of nearly 10 from peak to outer radius, although in some 
cases the drop of temperature from the peak value to that 
at r200 is consistent with the one reported in this paper. 

4.2 Modelling of the Chandra data out to rsoo 

The steepening of the radial profile beyond 400" is driven 
by the temperature of the last datapoint beyond r2oo- We 
also model the surface brightness and temperature pro- 
files of the Chandra data out to only 330", or a p proxi - 
mately r 5 oo, and find the best-fit IVikhlinin" et all (20061 ) 
model for the temperature profile reported in Figure [5] and 
Table [3] We measure a gas mass fraction of / g as(rsoo) = 
0.138 ± 0.006; if we add the mean stellar fraction as mea- 
su red by either iGiodini et all (|2009l ) (/* = 0.019 ± 0.002) or 
by I Gonzalez et aL | (l2007r> (/* ~ 0.012) assuming M(rsoo) = 



7.1 x 10 Mq, we find that Abell 1835 has an average baryon 
content within rsoo that is consiste nt with the cosmic ab un- 
dance of fl b /d M = 0.167 ± 0.007 (|Komatsu et al.ll201ll) at 
the 2-<7 level. As is the case in most clusters, especially re- 
laxed ones, the radial dist ribution of the gas m ass fraction 
increases with radius (e.g.. IVikhlinin et al.l l2006). 

We use this modelling of the data to measure rsoo, and 
to provide estimates for r2oo and the virial radius. The ex- 
trapolation of this model to 600" now falls above the mea- 
sured temperature profile, and the mass profile using hy- 
drostatic equilibrium is monotonic. This best-fit model is 
marginally compatible with the assumption of hydrostatic 
equilibrium. In fact, Table [4] shows that the extrapolated 



mass profile flattens around r2oo, with virtually no addi- 
tional mass being necessary beyond this radius to sustain 
the hot gas in hydrostatic equilibrium. Moreover, between 
rsoo and r2oo, all of the gravitational mass is accounted by 
the hot gas mass, i.e., no dark matter is required beyond 
rsoo- This extrapolation of the ^ rsoo data to the virial ra- 
dius therefore leads to a dark matter halo that is much more 
concentrated than the hot gas. 



5 ENTROPY PROFILE AND CONVECTIVE 
INSTABILITY AT LARGE RADII 

The Schwarzschild criterion for the onset of convective in- 
stability is given by the condition of buoyancy of an in- 
finitesimal blob of gas that is displaced by an amount dr, 
dpbiob < dp, where puob is the density of the displaced blob, 
assumed to attain pressure equilibrium with the surround- 
ing, and p is the density of ambient medium. If the blob is 
displaced adiabatically, using pressure P and entropy s as 
the independent thermodynamic variables in the derivatives 
of puob and p, the buoyancy condition gives 



dp 
ds 



ds>0 



(4) 



as condition for convective instability, i.e., a blob that is 
displaced radially outward will find itself in a medium of 
higher density and continue to rise to larger radii. Since 



(dp/ds) 



~p (dT/8P) s < (material is heated upon 



adiabatic compression), Equation [4] simply reads that a ra- 
dially decreasing entropy profile is convective unstable. 
An ideal gas has an entropy of 



S = vR ( T^mT- 



lnp + C 



(5) 



where v is the number of moles, R is the gas constant, and C 
is a constant. In astrophysi cal applications, it is customary 
(e.g. ICavagnolo et al.ll2009l ) to use a definition of entropy 
that is related to the thermodynamic entropy by an opera- 
tion of exponential and a constant offset, 



S = 



kT 

~ 273' 



(6) 



The entropy S defined by Equation [6] has units of keV 
cm 2 , and it is required to be radially increasing to main- 
tain convective equilibrium. Numerical simulations indi- 
cate that entropy outside the core is pred icted to increase 
with radius approxima tely as r ' or r ' (|Voit et al. I l2005l ; 
iTozzi fc Normanll200ll ). In Figure[l0]we show the radial pro- 
file of the entropy out to the outer radius of 10 arcmin, with 
a significant decrease at large radii that indicates an incom- 
patibility of the best-fit model with convective equilibrium. 
For comparison, we also show the entropy profile measured 
using the modelling of the data out to only rsoo, as described 
in Sec. 14.21 This entropy profile uses the shallower tempera- 
ture profile of Figure [9] and its extrapolation to larger radii 
remains non-decreasing, i.e., marginally consistent with con- 
vective equilibrium. 

The Schwarzschild criterion does not apply in the pres- 
ence of a magnetic field. For typical values of the thermody- 
namic quantities of the ICM, the electron and ion gyroradii 
are several orders of magnitude smaller than the mean free 
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Figure 7. Left: Best-fit Vikhlinin model for the projected temperature profile out to 600", with 90% confidence intervals. Right: Best-fit 
Vikhlinin model to the 0.7-2 keV surface brightness (model+background) profile. Emission beyond 600" is statistically consistent with 
the background, in blue is the extrapolation out to 1100". Prior removal of the stowed background caused the lower backgroud level in 
Figure [2] 



Table 3. Best-fit parameters for the Vikhlinin model using Chandra data out to 330" 



™e0 

(10" 2 cm- 3 ) 


(arcsec) 


P 


r s e n e0 2 
(arcsec) 


7 


a 


xL(d-o.f.) 


9.602±Bif| 


k 74 Q, + 0.373 


0.498±g : o°I 


Using Chandra data out to 330" 
119.8±J|;| 1.226±g;g|f 0.0 


3.0 


0.0 




9.763±8;ttS 


6.346±g;li 


0.488±g : o°l 


Using Chandra data out to 600" 
96.44±S:g 1.067±8;g^ 0.0 


3.0 


0.0 




To 
(keV) 


T 

(keV) 


fcool 

(arcsec) 


a coo l n at 
(arcsec) 


bt 


ct 




38.25±j^23 


3.0 


92.48±4o;i 


Using Chandra data out to 330" 

1.0 257.5±^ 3 ^ 0.0 


1.024±8;i| 


2.0 


39.0 (83) 


1 n 1 7+0.85 


3.0 


11.82±|;|J 


Using Chandra data out to 600" 
1.924±g;|g| 600.0 0.0 


2.800±°;|^ 


10.0 


106.4 (154) 



path for Coulomb collisions (e.g. ISarazi ni l 19881 ). even for a 
magnetic field of order 1 /xG, and the refore diffusion takes 
place primarily along field lines (e.g. IChandran fc Raseral 
l2007f ). There is strong evidence of magneti c fields in the 
central regions of clusters (e.g., radio halos, IVenturi et alj 
2008; Cas sano et alj [2006J ) . though it is not clear whether 
magnetic fields are u biquitous near the virial radius, as in 
the case of Abel l 3376 jBagchi et all20 06h In the presence of 
magnetic fields, IChandran fc Raseral (|2007l ) has shown that 
the condition for convective instability is simply dT/dR < 0. 

The Chandra data out to the virial radius therefore in- 
dicate that the ICM is convectively unstable, regardless of 
the presence of a magnetic field. In fact, in the absence of 
magnetic fields near the virial radius, Figure [TOl shows that 
Abell 1835 fails the standard Schwarzschild criterion, i.e., 
the entropy decreases with radius; in the presence of mag- 
netic fields, the negative gradient in the temperature profile 
alone is sufficient for the onset of convective instability (e.g., 
as discussed by IChandran fc Raseral 120071 '). Convective in- 
stabilities would carry hotter gas from the inner regions to- 
wards the ou ter region with in a few sound crossing times. 
As shown by ISarazinl (| 19881 ) . the sound crossing time for a 
10 keV gas is ~ 0.7 Gyr for a 1 Mpc distance, and an un- 



stable temperature gradient such as that of Figure [7] would 
be flattened by convection within a few Gyrs. Convection 
could in principle also result in an additional pressure gra- 
dient due to the flow of hot plasma to large radii, which can 
in turn help support the gas against gravitational forces. 



6 DISCUSSION AND INTERPRETATION 

In this paper we have reported the detection of X-ray emis- 
sion in Abell 1835 with Chandra that extends out to ap- 
proximately the cluster's virial radius. The emission can be 
explained by the presence of a cooler phase of the plasma 
that is dominant at large radii, possibly linked to the infall 
of gas from large-scale filamentary structures. We also inves- 
tigate the effects of clumping of the gas at large radii, and 
conclude that in principle a radial gradient in the clumping 
factor of the hot ICM can explain the apparent flattening of 
the entropy profile and the turn-over of the mass profile. 
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Table 4. Masses Calculated using Chandra data out to 330" , and Extrapolated out to noo 



A r A M gaa M total f gaa 

(arcscc) X 10 13 M X 10 14 M 

2500 164.9±|^ 

500 326. 6±il 

200 453.3±J|? 

100 570.9±2| : Si 



4.70±g;Jf 5.03±S|I 0.093±; 

10.75±g'^ 7.80±g : || 0.138±; 

15.36±g'|| 8.35±g : |f 0.184±; 

19.53±Sj : | 4 8.34±H2 0.234±; 




6.1 Detection of X-ray emission out to the virial 
radius 

The detection of X-ray emission out to a radial distance 
of 10 arcmin, or approximately 2.4 Mpc, indicates the 
presence of diffuse gas out to the cluster's virial radius. 
This is the first detection of gas out to the virial ra- 
dius with Chandra, matching other det ections obtained 
with Suzaku for ne a rby cl uste rs (e.g. Akamats u et alj 
20111; IWalker et all l2012bllal; ISimionescu et all l201ll; 



Burns et al.l 201p|; jKawaharada et alT 2010l : Bautz et alj 
20091 ; iGeorge et al l f2009l T Despite its higher background, 
Chandra provides a superior angular resolution to image 
and remove emission from unrelated sources. As can be 
seen from Figure [T] there are approximately 100 point-like 
sources that were automatically detected and removed, and 
we were also able to identify two low-mass clusters that 
are likely associated with Abell 1835. Chandra therefore 
has the ability to constrain the emission of clusters to the 
virial radius, especially for higher-redshift cool-core clusters 
for which the Suzaku point-spread function would cause 
significant contamination from the central signal to large 
radii. 

It is not easy to interpret the emission at the outskirts as 
an extension of the hot gas detected at radii ^ rsoo. In fact, 
as shown in Section |4. II the steepening of the temperature 
profile is incompatible with the assumption of hydrostatic 
equilibrium at large radii. We also showed in Section [5] that 
the gas has a negative entropy gradient beyond this radius, 
rendering it convectively unstable. Therefore, if the temper- 
ature profile of Figure[7]originates from a single phase of the 
ICM, convection would transport hotter gas towards the out- 
skirts, flattening the temperature profile within a few Gyrs. 



Cooling of the gas by thermal radiation cannot be responsi- 
ble for off-setting the heating by convection, since the cool- 
ing time (t coo i ~ kT 1 ^ 2 ^ 1 ) is longer at the outskirts than 
in the peak-temperature regions due to the higher density. 



6.2 Warm-hot gas towards the cluster outskirts 

A possible interpretation for the detection of emission near 
the virial radius and its steep temperature profile is the pres- 
ence of a separate phase at the cluster outskirts that is not in 
hydrostatic equilibrum with the cluster's potential. In this 
case, cooler gas may be the result of infall from filamentary 
structures that feed gas into the cluster, and the temper- 
ature of this warm-hot gas may in fact be lower than that 
shown in Figure[7](i.e., kT ~ 1.25 keV for the region ^ 450") 
if this gas lies in projection against the tail end of the hotter 
ICM. 

We estimate the mass of this putative warm-hot gas 
assuming that all of the emission from the outermost re- 
gion is from a uniform density gas seen in projection. This 
assumption may result in an overestimate of the emission 
measure; in fact, the extrapolation of the gas density profile 
in the hydrostatic or convective scenarios may yield a sig- 
nificant amount of emission in the last radial bin. We were 
unable to perform a self-consistent modelling of the emis- 
sion in the full radial range, since the low signal-to-noise 
ratio does not allow a two-phase modelling in the last radial 
bin. In this simple uniform density warm-hot gas scenario, 
the gas is in a filamentary structure of length L and area 
A = Tv{R 2 out - R 2 n ), where R out = 600" and R in = 450"; 



this is the same model also considered in lBonamente et al.l 
(2005) for the cluster Abell SI 101. Since the length L of the 
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Figure 9. Temperature and gas mass fraction profiles measured from a fit to the Chandra data out to 330", and extrapolation of the 
best-fit model out to 600" . 
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Figure 10. Deprojected entropy profiles using the full Chandra data out to 600" (left, see Scction l4.ll . and using only data out to rsoo 
(right, see Section 1 1.211 . 



filament along the sightline is unknown, we must either as- 
sume L or the electron density n e , and estimate the mass 
implied by the detected emission. The emission integral for 
this region is proportional to 

K = ^mrw n ^ (7) 

where K is measured in XSPEC from a fit to the spectrum, 
Da is the angular distance in cm, z is the cluster redshift, 
and the volume is V = A x L. For this estimate we assume 
for simplicity that the mean atomic weights of hydrogen 
and of the electrons are the same, /J, e = Using the best- 
fit spectral model with kT = 1.26 ± 0.16 keV, we measure 
K = 1.05 ±0.13 x 10 -4 . If we assume a filament of length 
L — 10 Mpc, then the average density is n e = 2.4±0.3 cm~ 3 , 
and the filament mass is 4.6 ±0.6 x 10 13 M e . Alternatively, 
a more diffuse filament gas of n e = 10~ 5 cm" 3 would require 
a filament of length L = 58 ± 8 Mpc, with a mass of 1.1 ± 
0.2 x 10 14 Mq, comparable to the entire hot gas mass within 
r2oo ■ The fact that a lower density gas yields a higher mass is 
given by the fact that, for a measured value of K we obtain 
n e oc L~ 1 ^ 2 , and therefore the mass is proportional to L 1 ^ 2 . 
For comparison, the gas mass for this shell inferred from the 



standard analysis, i.e., assuming that the gas is in the shell 
itself, is ~ 3 x 10 13 M Q , as can be also seen from Table [4] 

If the gas is cooler, then the mass budget would increase 
further. In fact, the bulk of the emission from cooler gas falls 
outside of the Chandra bandpass, and for a fixed number of 
detected counts the required emission integral increases. We 
illustrate this situation by fitting the annulus to an emission 
model with a fixed value of KT = 0.5 keV, which result in 
a value of K = 1.88 ± 0.24 x 10" 4 (the fit is significantly 
poorer, with A^ 2 = ±10 for one fewer degree of freedom). 
Accordingly, the filament mass estimates would be increased 
approximately by a factor of two. 

A warm-hot phase at T ^ 10 7 K is expected to 
be a significant reservoir of baryons i n the universe (e.g. 
ICen fc Ostrikerlll999l ; iDave et al.(l200ll ). Using the ROSAT 
soft X-ray Position Sensitive Proportional Counter (PSPC) 
detector - better suited to detect the emission from sub- 
keV plasma - we have already reported the detection of 
a large-scale halo of emission around the Coma cluster 
out to ~ 5 Mpc. well beyon d the cluster's virial radius 
l|Bonamente et all 120031 , [2009). It is possible to speculate 
that the high mass of Abell 1835, one of the most luminous 
and massive clusters on the Bright Cluster Survey sample 
|Ebeling et a.1.1 Fl998) . is responsible for the heating of the 
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infalling gas to temperatures that makes it detectable by 
Chandra, and that other massive clusters may therefore pro- 
vide evidence of emission to the virial radius with the Chan- 
dra ACIS detectors. The infall scenario is su pported by the 
Herschel observations of iPereira et al.l (|2010l ). who measure 
a galaxy velocity distribution for Abell 1835 that does not 
appear to decline at large radii as in most of the other clus- 
ters in their sample. A possible interpretation for their data 
is the presence of a surrounding filamentary structure that 
is infalling into the cluster. 

6.3 Effects of gas clumping at large radii 

Masses and entropy measured in this paper assume that 
the gas has a uniform density at each radius. To quantify 
the effect of departures from uniform density, we define the 
clumping factor C as the ratio of density averages over a 
large region, 



C = 



(nl) 

K> 2 



(8) 



with C ^ 1. Clumped gas emits more efficiently than gas of 
uniform density, and the same surface brightness / results 
in a lower estimate for the gas density and mass, 



I oc / < n» > dl — 



< n, 



> 2 Cdl, 



(9) 



where I is a distance along the sightline. From Figure [10] 
we see that the entropy drop from approximately 400" 
to 600" would be offset by a decrease in n^ 3 by a fac- 
tor of 3, or a decrease in n e by a factor of 5. We there- 
fore suggest that a clumping factor of C ~ 25 at 600" 
would in principle be able to provide a flat entropy profile, 
and even higher clumping factors would provide an increas- 
ing entropy profile in better agreemen t with theory (e.g. 
IVoit et al.ll2005l; iTozzi fc Normanll200ll ). Numerical simula- 
tions by iNagai fc Laul (|201ll ) suggest values of the clump- 
ing factor C < 3 near r2oo, with signifi cantly higher c lump - 
ing possible at larger radii. Use of the INagai fe Laul (120111 ) 
model in the analys is of a large sample of galaxy clusters by 
lEckert et al. (2012) results in better agreement of observa- 
tions with numerical simulations. 

Clumping can also affect the measurement of hydro- 
static masses. In particular, gas with an increasing radial 
profile of the clumping factor could result in a steeper gra- 
dient of the density profile, when compared with what is 
measured assuming a uniform density. According to Equa- 
tion[2] this would result in larger estimates of the hydrostatic 
mass, in principle able to reduce or entirely offset the ap- 
parent decreas of M(r) reported in Figure [8] We therefore 
conclude that a radial increase in the clumping of the gas 
can in principle account for the apparent decrease of the 
mass profile and of the entropy profile reported in this pa- 
per (Figures [8l and [10)) . and therefore it is a viable scenario 
to interpret our Chandra observations. Clumping of the gas 
at large radii h as also been suggested based on Suzaku ob- 
servations (e.g.. Isimionescu et al.llioill ). 



7 CONCLUSIONS 

In this paper we have reported the detection of emission 
from Abell 1835 with Chandra out to the cluster's virial ra- 



dius. The cluster's surface brightness is significantly above 
the background level out to a radius of approximately 10 
arcminutes, which correspond to ~2.4 Mpc at the cluster's 
redshift. We have investigated several sources of systematic 
errors in the background subtraction process, and deter- 
mined that the significance of the detection in the outer 
region (450-600") is ^ 4.7 a, and the emission cannot be ex- 
plained by fluctuations in the background. Detection out to 
the virial radius is also impli ed by the XMM-New ton tem- 
perature profile reported bv lSnowden et all (|2008l ) . 

The Chandra superior angular resolution made it 
straightforward to identify and subtract sources of X-ray 
emission that are unrelated to the cluster. In addition to a 
large number of point sources, we have identified X-ray emis- 
sion from two low-mass clusters that were sel ected from the 
SPSS data, MAXBCG J210.31728+0 2. 75364 (iKoester et all 
l2007f) and WHL J140031.8+025443 (|Wen et al.ll2009l ). The 
two clusters have photometric and spectroscopic redshifts 
that make them likely associated with Abell 1835. These are 
the only two SDSS-selected clusters that are in the vicinity 
of Abell 1835. 

The outer regions of the Abell 1835 cluster have a sharp 
drop in the temperature profile, a factor of about ten from 
the peak temperature. The sharp drop in temperature im- 
plies that the hot gas cannot be in hydrostatic equilibrium, 
and that the hot gas would be convectively unstable. A pos- 
sible scenario to explain the observations is the presence of 
warm-hot gas near the virial radius that is not in hydro- 
static equilibrium with the cluster's potential, and with a 
mass budget comparable to that of the entire ICM. The data 
are also consistent with an alternative scenario in which a 
significant clumping of the gas at large radii is responsible 
for the apparent negative gradients of the mass and entropy 
profiles at large radii. 
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